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We briefly review the theory of Bose-Einstein condensation in the two-dimensional 
trapped Bose gas and, in particular the relationship to the theory of the homogeneous 
two-dimensional gas and the Berezinskii-Kosterlitz-Thouless phase. We obtain a phase 
diagram for the trapped two-dimensional gas, finding a critical temperature above which 
the free energy of a state with a pair of vortices of opposite circulation is lower than 
that for a vortex-free Bose-Einstein condensed ground state. We identify three distinct 
phases which are, in order of increasing temperature, a phase coherent Bose-Einstein 
condensate, a vortex pair plasma with fluctuating condensate phase and a thermal Bose 
gas. The thermal activation of vortex-antivortex pair formation is confirmed using finite- 
temperature classical field simulations. 
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1. Introduction 

It is now over a decade since Bose-Einstein condensation (BEC) in a dilute atomic 
gas was first realised.^ In this time remarkable progress has been made in the con- 
trol and manipulation of these ultra-cold gases. One possibility offered, through the 
use of light fields in the form of the optical lattice,^ by this control is to reduce the 
dimensionality of the system. By increasing the confining potentials in one or two 
directions it is possible to freeze out these degrees of freedom creating, effectively, 
two-^ or one-dimensional^ gases. The two-dimensional (2D) gas is of particular in- 
terest. In general 2D systems often display unique properties, often remarkably 
different from those of the comparable bulk system. For example, the 2D electron 
gas, subjected to a perpendicular magnetic field, displays the rich phenomena of 
the, essentially single particle, quantum Hall effect and, the even more remark- 
able, many body fractional quantum Hall effect^ with the physical manifestation of, 
the uniquely two-dimensional, fractional statistics^ in the quasiparticle excitations. 
More specifically, it is well known that the uniform 2D Bose gas does not un- 
dergo BEC^'^ at finite-temperature. The interacting homogeneous 2D gas does still 
undergo a normal-superfluid phase transition^'^"'^^ at finite-temperature however, 
with an order parameter which is only locally phase coherent - the Berezinskii- 
Kosterlitz-Thouless (BKT) phase. Experimental evidence for this phase transition 
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has been demonstrated in liquid helium thin film^^, superconducting Josephson- 
junction arrays^^ and in spin-polarized atomic hydrogen^^. 

In the trapped 2D Bosc gas the density of states is modified sufficiently 
that, for the ideal gas, condensation into a single state can occur at finite 
temperature.^^ It has been demonstrated however that long- wavelength phase fluc- 
tuations may destroy the global phase coherence of the condensate at a tem- 
perature still below the ideal gas critical temperature, To for BEC with only a 
quasi-condensate^^'^'^'^^'^^''^^''^^''^'^ present. The connection between such a quasi- 
condensate and the correlated vortex-pair plasma of the BKT phase proves unclear 
however and forms the basis of the material we attempt to address in the remainder 
of this paper. 

2. Free Energy of the Vortex- Antivortex Pair 

In order to determine at what temperature thermal activation of vortex-antivortex 
pairs becomes thermodynamically favourable, we consider the Helmholtz free en- 
ergy, F = E — TS of bosonic particles of mass m spatially confined in 2D and 
via the harmonic potential ^(r) = mco'^r'^ /2, where u)± is the radial trap frequency. 
We assume the existence of a macroscopic ground state condensate wavefunction 
ip{r), which also serves to describe the order parameter of the system in the BEC 
phase. Isolating the energy contribution due to the condensate, Eq{T), the total in- 
ternal energy of the system may be written as E{T) = Eo{T) + E{T). The required 
condensate energy Eo{T) is determined by the functional 

= / (I^IVV-Wr + V(r)|V(r)r + f l^(r)l^) dr (1) 

where g is the constant coupling parameter for the particle interactions, and 
Nq = J |'0(r)pdr denotes the number of particles in the condensate. The internal en- 
ergies for topologically distinct order parameter configurations are then calculated 
computationally.^^ The entropy contribution to the free-energy difference is due to 
the multiplicity of order parameter configurations containing a pair of vortices. The 
statistical weight W = , is obtained by allowing one vortex to reside any- 

where within the Thomas-Fermi radius, the partner vortex then having 27r available 
nearest neighbour sites. Assuming the radius of the area occupied by each vortex 
to be of the order of the healing length, ^ = ■\/fi^/2m/i, we can obtain an expres- 
sion for the configurational energy difference AS{T) = ks In (STr/x^/fi^w^^). The free 
energy for each configuration can now be evaluated, yielding a critical temperature 
for the thermal activation of vortex-antivortex pair creation. For the experimental 
parameters of the Oxford experiments^ we obtain a critical temperature Tc ~ O.STo 
although the precise value is dependent upon the total number of atoms and de- 
creases monotonically with A''. Of course taking the thermodynamic limit in this 
case, which one would need to do if one were to try to claim that this is a phase 
transition, is slightly complicated by the presence of the trap. 
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This indicates that there are three distinct regions for the harmonically trapped, 
ultra-cold, dilute Bose gas. At low temperatures there is a phase coherent BEC 
with a coherence length comparable with the system size as shown by Gies and 
Hutchinson.-^^ At higher temperatures, where the coherence length reduces and off- 
diagonal order decays algebraically, it becomes thermodynamically favourable to 
form correlated vortex-antivortex pairs, which we take as a signature of a BKT-like 
phase. At temperatures approaching the ideal gas critical temperature, Tq the off- 
diagonal order decays exponentially. At this point the vortex pairs unbind (although 
this is not present in our free energy model) and superfluidity is destroyed yielding 
a thermal gas. This is what is traditionally referred to as the BKT transition. 



3. Classical Field Simulations 

To confirm these conclusions finite temperature classical field simulations^®'^^ were 
performed. In the classical field approach the system is divided into classical and 
incoherent regions, determined by the occupation of single particle modes. Highly 
occupied (classical) modes are described by the projected Gross-Pitaevskii equation 



ih 



dt 



1 



*(r,t)+!,-p{Wr,t)lV(r,t)} (2) 



where the projector, restricts evolution of the classical field to within its subspace. 
The low occupation, incoherent modes are described using the scmiclassical Hartrcc- 
Fock approximation, with the classical and incoherent regions being taken to be in 
thermal equilibrium with one another. At the lowest temperatures no vortices are 
present in the system. These first emerge, as the temperature is increased, in the 
low density regions at the edge of the cloud. At higher temperatures it becomes 
possible for vortex-antivortex pairs to nucleate closer to the trap centre. This is 
demonstrated in Fig.(l) which is at a temperature of O.86T0. Panel (a) shows a 
single frame of the full simulation with vortices and antivortices shown by the 
black crosses and white dashes. These are identified through the singularity in the 
phase as determined from the simulation. The region containing a slight density dip 
(due to thermal fluctuations) bounded by the box is then blown up into panel (b). 
Subsequent time frames are then shown for this expanded region in panels (c), (d) 
and (e). It can clearly be seen how a vortex-antivortex pair forms in the low density 
fluctuating region, then separates and stabilises. At later time this pair goes on to 
annihilate again. The temperature at which vortex-antivortex pairs begin to form in 
the classical field simulations is entirely consistent with the free energy calculations 
described above and a strong corroboration of the interpretation of the regime with 
only algebraic off-diagonal order as a BKT-like phase. Further details of the classical 
field simulations for the 2D gas and comparison with recent experiments^^ can be 
found in Ref. 29. 
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Fig. 1. Panel (a) shows a single frame of the full simulation with vortices and antivortices shown 
by the black crosses and white dashes. The region bounded by the box is blown up into panel 
(b). Subsequent time frames are shown for this expanded region in panels (c), (d) and (e) clearly 
showing the vortex-antivortex pair nucleation. 

4. Conclusions 

In conclusion we have demonstrated, both using a free energy argument and us- 
ing classical field simulations that there exist three distinct regions in the phase 
diagram of the ultra-cold, harmonically confined 2D Bose gas. At very low tem- 
peratures compared to the ideal gas condensation temperature Tq the gas forms a 
phase coherent BEC. At temperatures of order 0.5To it becomes thermodynami- 
cally favourable for correlated vortex-antivortex pairs to form. The system remains 
superfluid with algebraic off-diagonal order. This corresponds to a BKT-like su- 
perfluid phase. At higher temperatures, approaching Tq, the vortex pairs unbind, 
the off-diagonal density matrix decays exponentially and the superfluid phase is 
destroyed yielding a thermal gas. This upper transition is what corresponds to the 
BKT transition in the uniform system. 
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